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ABSTRACT 


Rotordynamical analyses are typically performed using rigid disk models. Stud- 
ies of rotor models in which the effects of disk flexibility have been included indicate 
that it may be an important effect for many systems. This work addresses this issue 
with respect to the Space Shuttle Main Engine high pressure turbo-pumps. Finite 
element analyses have been performed for a simplified free-free flexible disk rotor 
model and the modes and frequencies compared to those of a rigid disk model. The 
simple model was then extended to a more sophisticated HPTOP rotor model and 
similar results were observed. Equations have been developed that are suitable for 
modifying the current rotordynamical analysis program to account for disk flexibil- 
ity- Some conclusions are drawn from the results of this work as to the importance 
of disk flexibility on the HPTOP rotordynamics and some recommendations are 
given for follow-up research in this area. 
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NOMENCLATURE 


h = disk thickness coordinate 
I a ,» = axial moment of inertia 
I t ,i = transverse moment of inertia 
= ///(h 2 rp)drd*/dh 
M = mass of rotor 
m = mass matrix of rotor 
< 7 , = i th generalized coordinate 
r = disk radial coordinate 

x iy Viy z % — Translational coordinates of rotor disk hub 
X = Lateral translational coordinate of rotor disk hubs 
& = i th generalized coordinate for non — rotating system 
$ = mode shape 
^ = mode shape 
A = mode shape 
v = disk angular coordinate 
p = material density 

r = W<AH*> 

A = {*} T {/„}{¥} 

ft = Rotor speed 

0 yy $ My 9 X — Euler angles 
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INTRODUCTION 


In the modeling of any physical device or process, certain assumptions and 
restrictions must be made. It is important to carefully assess their validity in order 
to determine the accuracy and range of applicability of the mathematical model. 
The current procedure for analyzing the rotordynamics of the Space Shuttle Main 
Engine Turbo- pumps consists of coupling the free-free rotor and housing modes 
with constraints to produce a model for the complete turbo-pump. This model is 
then used in stability analyses and for time response simulations. 

The standard practice in obtaining the free-free rotor modes is to neglect the 
effects of rotor disk flexibility. Research into this area has indicated that disk 
flexibility may play an important role in the rotordynamical behavior of turbo- 
machinery. If rotor disk flexibility has a significant effect on the rotordynamics of 
the SSME turbo-pumps, significant errors could be introduced into analyses. So, it 
is important to assess such effects and develop means of accounting for it in analysis 
procedures. This work attempts to address this issue by examining some simplified 
finite element rotor motors and developing analytical methods of dealing with disk 
flexibility effects. 
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OBJECTIVES 


(1.) Develop finite element rotor models to evaluate the influence of disk flexibility 
on the rotor free-free modes. 


(2.) If justified by (1.), develop methods of modifying the current rotordynamical 
analysis program to account for rotor disk flexibility. 


3.) Develop recommendations for follow-up research. 
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PREVIOUS RESEARCH 


There is a large body of research that is available on the dynamical behav- 
ior of rotating flexible bodies. A classic work by Lamb and Southwell presented 
a discussion of the vibrational behavior of a spinning disk. 1 This work served as 
the basis for much further work in this area. Likins performed a study of math- 
ematical modeling of spinning elastic bodies for use in modal analyses. 2 Wilgens 
and Schlack investigated the dynamical behavior of a flexible beam attached to a 
rotating snaft. 3 Brown and Schlack extended this work to a study of the stability 
of a spinning body. 4 Pringle presented a method for examining the dynamical be- 
havior of a system with connected moving parts. 6 Laurenson discussed methods for 
performing modal analysis on rotating flexible structures. 6 Meirovitch presented an 
analytical study of discretization methods for flexible gyroscopic systems and drew 
conclusions concerning the appropriateness of various techniques. 7 

There is also a fairly large body of work documented in the literature concern- 
ing studies of disk flexibility on rotors and turbomachinery. Vance studied several 
rotor systems, performing a combined experimental - analytical investigation. 8 His 
analytical model was configured so that the effects of disk flexibility could be ac- 
counted for in an approximate fashion. These studies indicated that, for most cases, 
the inclusion of rotor disk flexibility could significantly improve the correlation be- 
tween experimentally measured rotor free-free natural frequencies and calculated 
values. Klompas developed a technique for studying shaft whirling which included 
the effects of flexible disks and blades. 0 Klompas extended his study further by 
including the effects of disk flexibility in a study of the unbalance response of a 
turbo-machine. 10 This work was primarily aimed at a study of the effects of blade 
loss. Palladine and Rosettos developed a finite element method for examining the 
effects of flexibility on the behavior of a rotor. 11 Wilgen and Schlack investigated 
the effects of disk flexibility on shaft whirl stability using Liapunov techniques. 12 
The resulting procedure is analytically very nice, but for complicated disk shapes 
and multi-disk systems the method could quickly become intractable. Dopkin and 
Shoup performed a study of the effects of disk flexibility on the resonant frequencies 
of an axisymetric rotating shaft. They found that the effects of disk flexibility may 
significantly reduce the rotor resonant speeds and that this effect was particularly 
pronounced at low rotor speeds. 13 Shahab and Thomas studied finite element mod- 
els of single and multi-disk rotor systems and compared the results to experimental 
models. 14 This study indicated that coupling effects between the shaft and disk 
modes can have a significant effect on the dynamical behavior of a rotor. Sakata, 
Aiba, and Ohnabe studied the transient vibration behavior of a rotor subjected to 
a blade loss and included the effects of disk flexibility. 16 
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FINITE ELEMENT ANALYSES 


The first objective of this research effort is aimed at evaluating the influence 
of disk flexibility on the rotor free-free modes. The most straight-forward way 
of studying this effect is by developing finite element rigid and flexible disk rotor 
models and comparing the rotor free-free modes. 

Examination of schematics of the HTOTP and HTFTP rotor disk configura- 
tions reveals that the second turbine stage of the HTOTP has the thinnest disk 
configuration and is probably most likely to exhibit flexibility effects. In order to 
establish what types of behavior might be expected for a rotor with a flexible disk, 
a simple model was first examined. A configuration was selected that is approxi- 
mately that of the HTOTP rotor shaft with a disk attached to represent the second 
turbine stage. This model consists of a 0.051 m. diameter, 0.59 m. steel shaft. 
A 0.239 m. diameter, 0.0161 m. thick disk was attached to the shaft with its 
center point 0.528 m. from the end of the shaft. For the flexible disk study, the 
properties of steel were used for the disk. For the rigid disk study, the density and 
poisson’s ratio of steel and a modulus of rigidity and modulus of elasticity three 
orders of magnitude above those of steel were used as the material properties. The 
free-natural frequencies for the two models are compared in Table 1 and the mode 
shapes are illustrated in Figure 1. 


Table 1 : Free-Free Natural Frequencies 
For Simple Rotor Models 


Natural Frequencies (Hz) 

Mode 

Flexible Disk 

Rigid Disk 

First Bending 

487 

502 

Second Bending 

1003.0, 1997.0 

1225.0 

Third Bending 

3393.0 

2453.0 


Note that the first rotor bending mode is not significantly affected by the disk 
flexibility. However, the second and third bending modes are strongly influenced. 

In order to better evaluate the influence of disk flexibility, it was decided to 
develop a second, more sophisticated rotor model of the HPTOP. A finite element 
code was developed based on work by Muller 10 . The model was developed using 
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ANSYS, a finite element analysis package. The types of beams and rigid masses 
used are the same as those used by Muller except the second turbine has been 
replaced by a 0.239 m. diameter, 0.0161 m. thick disk, and rigid inertias I xx = 
5.28x10 3 kg - m 2 , I yy = 2.96rl0~ 3 kg — m 2 , and I„ = 2.96xl0” 3 kg - m 2 . As in 
the previous study, the material properties of steel are used for the flexible disk case. 
For the rigid disk, the density and poisson’s ratio of steel are used with a modulus 
of elasticity and modulus of rigidity increased by three orders of magnitude. The 
resulting model was examined using an eigenanalysis to determine the free-free 
rotor modes. 


Table 2 : Free- Free Natural Frequencies 
For HPOTP Rotor Models 


( Natural Frequencies (Hz) 

Mode 

Flexible Disk 

Rigid Disk 

First Bending 

462.2 

467.6 

First Torsional 

913.8 

944.3 

Second Bending 

926.7, 1652.0 

1072.0 

Second Torsional 

1255.2 

1274.0 

Third Bending 

1773.6 

1983.0 

Third Torsional 

2636.7 

2486.0 


For the first bending modes, the mode shapes and frequencies for the two cases 
match closely. Similar behavior is observed for the first torsional modes. For the 
second bending modes, the rigid disk system has a single mode shape. The flexible 
disk model exhibits two second bending modes. The first is characterized by motion 
of the disk in-phase with the hub. The second is characterized by motion of the 
disk out-of-phase with the hub. The in-phase motion serves to effectively reduce 
the second bending mode natural frequency and the out-of-phase motion serves to 
effectively increase it. 
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ANALYTICAL STUDY 


The results of the finite element studies indicate that rotor disk flexibility 
can significantly alter the free-free rotor modes and frequencies. As a result, it 
is appropriate to develop techniques for accounting for such effects in the current 
rotordynamical analysis program so that the significance of disk flexibility can be 
evaluated for the complete turbo-pump model. 

The equations of motion for two formulations of a flexible rotor - flexible disk 
model are presented. Each of the approaches presented in this study assume that 
the deflection of the rotor disk is primarily in the lateral direction. The beam is 
assumed to be axially and torsionally rigid, which implies that the disk hubs move 
together as a rigid body in the lateral direction. In addition, it is assumed that the 
beam is axially stiff so that the axial hub position of i ih disk is = X y for all i. 
For each of the following developments, the rotor is considered to consist of a series 
of flexible disks. The equations are formulated using a Lagrangian formulation. 


Rotation Sequence 


9 y i about y 
0 Mi about zt 
0 x i about x 

v about i (disk angular coordinate) 


Formulation Using Flexible Rotor - Flexible Disk Modes 

Position of Arbitrary Point 
In Inertial Coordinate System 

First, obtain the position vector P{ for an arbitrary point on the i* h rotor disk. 
Pi = P x ,iX + P yi y + P ti z 


where 



Py.i = Yi + hsin{9 Zti ) + rcos(i/ + 0 Xi< )cos(0, (< ) 

P z i — Z, + rsin(^ + 9 Xti ) cos(9 Vii ) - /tco8(0 S|< )sin(0 v ,{) 

+rcos(i/ + 9 X i) sin(0 y< ) sin(0,,i) 

Velocity of Arbitrary Point 

Differentiation of the position vector with respect to time yields a velocity 
vector for the arbitrary point. 

?i=V a , i i + V 9 , i y + V Mii i 

where 

V,,i = ti, 

V y i = Y~ 9 z i hcos(9 z i ) -rnsin(i/ + ff Ii j)cos(tf, i< ) —r9 z>i cos(i/ + 0*,<)sin(0.,,) 

V, ,- = Zi + rfJ cos(i/ + 9 a<i ) cos(0 y> j) - r0 y> , sin(v + 0*,, ) sin(0 y> , ) 

+h9 z i sin(0*.,) sin(0 yi< ) - h9 v>i cosftf,,,) cos(0 v>i ) 

-rflsin(i/ + 0 x ,,)sin(# Vl <)8in(0,,t) +r0 Vi< cos(i/ + tf« i< )cos(tf Vi , )sin(tf li j) 
+r9 Z 'i cos(i/ + 0 Ii ,)8in(tf y>< )cos(^, i< ) 

Kinetic Energy 

Express the kinetic energy of each disk as 
Ti = \Sf f(V*S + V*/ + V M S)prdrdvdh 

The total kinetic energy of the rotor is then 
T = E< 7; 

T =\ f f j {{u} T {<i})rdrdudh 
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+i{*.} T [ih](«.} +i{<U r H,]{«„> + jn{#,nu]{#.} + jO{#.) T M{*,} 
+jn J M +|{Y) T H{r } 1 +i(z} T (m]{i | 1 -in*{«.} T [i.]{»,} 

where /*,, = f f f (h 2 pr)drdvdh 


{4} = 


h,l 

0 

0 

0 

ha 

0 

0 

0 

0 



(Z) 




\u n J 


For a non-spinning rotor, the position of a point on the disk can be expressed 
in terms of the free-free rotor disk modes. 


{u} = {h} + {X}+^ 



where X is the hub axial position. 

For a spinning rotor, it is necessary to transform this relation to account for 
the rotor spin. Notice that u* is defined in terms of the inertial reference frame. In 
order to express in terms of the free-free disk modes, one can make the following 
coordinate transformation. 

(tv) = cos(fi£){g„} + sin(fi/){g,} 
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{£,} = cos(nt){g,} - sin(nf){ 9v } 


<»> = {*} + <*}+(*' 


o W cos(nt) sin(n^) \ ( qy\ 
A v ) \ - sin(fit) cos(fi*) ) \q z ) 


Notice that since the rotor is axisymmetric, 

/ f f{[ A y] T {^y} r P) dr(ivdh = f J f{[^z] T [^M]rp)drdi/dh 

Let [A] - J f J ({A y ] T [A v ]rp)drdi/dh 

and 

t = i{«») r [A|«») +^«.} t [a]{?.> -jn«.Hr|{ fc } -tn{j„} r |r]{ ? ,} 

+i[n|[i.![n| +iMX> +i{«„> r {*} T {m){*)a) 

+ ^{q y } T [A]{q.} -~n{q,} T [A]{q,} + ^n{q,} T [A]{q,} 
4!){ q ,r[A]<q,> +~H a {?,} r |A]{?„} +ilJ J {».} T |A|{,.} 
+j{fc> r |A|<fc>+|M.} r |A]«.} 

r = {*> T (/„]{*} 

a = WWW 


Potential Energy 
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5 10 


ft 



0 




n2 


n 


= (?» 



V o 

0 . 

• • 

[V] 

Ill 

Of 

0 

Ml. 

,(nt) 

- sin(nt) \ / 

'Ml 

0 ' 

(nt) 

cos(nt) J \ 

< 0 

ML 


Ill 

, 0 

o > 

ML 


<3- 


l(fii) 


ain(nf) \ / g„ \ 
cos(nt ) ) \q z ) 


where [w*] represents the squared natural frequencies of the flexible rotor - 
flexible disk. 


Equations of Motion 


{$} r {m}{$}^ + [Aj{g,} +[A }{q v } +2n[A]{q.} +In 3 [r]{q y } -fi 2 [A }{q v } + 

( w nK?y} = {0} 


{$} T {m}{$}g; + [A] {q M > +[A]{g;} -2n[A]{q y } +in 2 [r]{q t } +n 2 [A]{g,} + 
]{?,} = {0} 


{X} = {0} 
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CONCLUSIONS AND RECOMMENDATIONS 


From the results of the finite element analyses, it is clear that rotor disk flex- 
ibility can significantly alter the rotor free-free modes and frequencies. While the 
first rotor bending mode is not strongly affected by disk flexibility, the second and 
third bending modes are significantly altered. In fact, two second bending modes 
are identified. The first is associated with in-phase motion of the disk with the 
rotor and the second is associated with out-of-phase motion of the disk with the 
rotor. 

Equations have been developed for accounting for disk flexibility in a rotor 
model. Particular emphasis has been placed on obtaining equations that are suitable 
for incorporation into the current used rotordynamical analysis program. 

From these results, the following conclusions and recommendations have been 
drawn. 


1. ) The rotordynamical analysis program should be modified to account for disk 

flexibility. 

2. ) The revised program should be tested with modal data from a simplified rotor 

finite element model. 

3. ) If results warrant, develop full-scale finite element models of the SSME turbo- 

pump rotors and use the resulting modal data in the revised rotordynamical 
analysis program. 

4. ) Compare these results to the responses predicted for a rigid disk rotor and 

evaluate the influence of rotor disk flexibility on the SSME turbo-pumps. 

5. ) In order to develop further physical insight into the effect of rotor disk flex- 

ibility, construct appropriately scaled rotor models and study their responses 
using a rotor test kit. 

6. ) Relate the results of these studies to actual observed behavior of the SSME 

turbo-pumps in order to gain physical insight and understanding of their ro- 
tordynamical behavior. Such understanding could enhance failure analysis . 
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FIGURE l.a: RIGID DISK MODEL - FIRST BENDING MODE 

FREQUENCY = 502 Hz 



FREQUENCY = 487 Hz 


FIGURE 1 : SIMPLE ROTOR MODELS - FIRST BENDING MODES 
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FIGURE 3. a: 


RIGID DISK MODEL - THIRD BENDING MODE 

FREQUENCY = 2453 Hz 



FIGURE 3.b: FLEXIBLE DISK MODEL - THIRD BENDING MODE, DISK OUT -OF -PHASE 

FREQUENCY = 3393 Hz 


FIGURE 3: SIMPLE ROTOR MODELS - THIRD BENDING MODES 


FIGURE 4. a: RIGID DISK MODEL - FIRST BENDING MODE 

FREQUENCY = 467.6 Hz 



FIGURE 4 . b : FLEXIBLE DISK MODEL - FIRST BENDING MODE 

FREQUENCY = 462.2 Hz 


FIGURE 4: HPOTP ROTOR MODELS - FIRST BENDING MODES 



FIGURE 5. a 


RIGID DISK MODEL - 
SECOND BENDING MODE 



IN-PHASE 

FREQUENCY = 926.7 Hz 



FIGURE 5.c 



FLEXIBLE DISK MODEL - 
SECOND BENDING MODE, DISK OUT-OF-PHASE 
FREQUENCY = 1652.0 Hz 


FIGURE 5: HPOTP 


ROTOR MODELS - SECOND BENDING MODES 



FIGURE 6. a: RIGID DISK MODEL - THIRD BENDING MODE 



FIGURE 6 . b : FLEXIBLE DISK MODEL - THIRD BENDING MODE 

FREQUENCY = 1773 Hr 


FIGURE 6: HPOTP ROTOR MODELS - THIRD BENDING MODES’ 


XI-19 


ORIGINAL PAGE IS 
OF POOR QUALITY 



ORIGINAL PAGE IS 
OF POOR QUALITY 



FIGURE 7: HPOTP FLEXIBLE DISK ROTOR MODEL - DISK MODE 
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